Point-Source Power in 3 Year Wilkinson Microwave Anisotropy Probe Data

Using a set of multifrequency cross spectra computed from the 3 year WMAP sky maps, we fit for the unresolved point-source contribution. For a white-noise power spectrum, we find a Q-band amplitude of A = 0.011 ± 0.001 μK^2 sr (antenna temperature), significantly smaller than the value of 0.017 ± 0.002 μK^2 sr used to correct the spectra in the WMAP release. Modifying the point-source correction in this way largely resolves the discrepancy that Eriksen et al. found between the WMAP V- and W-band power spectra. Correcting the co-added WMAP spectrum for both the low-l power excess due to a suboptimal likelihood approximation—also reported by Eriksen et al.—and the high-l power deficit due to oversubtracted point sources—presented in this Letter—we find that the net effect in terms of cosmological parameters is an ~0.7 σ shift in n_s to larger values. For the combination of WMAP, BOOMERANG, and ACBAR data, we find ns = 0.969 ± 0.016, lowering the significance of n_s ≠ 1 from ~2.7 σ to ~2.0 σ

Measuring Planck beams with planets

Aims. Accurate measurement of the cosmic microwave background (CMB) anisotropy requires precise knowledge of the instrument beam. We explore how well the Planck beams will be determined from observations of planets, developing techniques that are also appropriate for other experiments. Methods. We simulate planet observations with a Planck-like scanning strategy, telescope beams, noise, and detector properties. Then we employ both parametric and non-parametric techniques, reconstructing beams directly from the time-ordered data. With a faithful parameterization of the beam shape, we can constrain certain detector properties, such as the time constants of the detectors, to high precision. Alternatively, we decompose the beam using an orthogonal basis. For both techniques, we characterize the errors in the beam reconstruction with Monte Carlo realizations. For a simplified scanning strategy, we study the impact on estimation of the CMB power spectrum. Finally, we explore the consequences for measuring cosmological parameters, focusing on the spectral index of primordial scalar perturbations, n_s. Results. The quality of the power spectrum measurement will be significantly influenced by the optical modeling of the telescope. In our most conservative case, using no information about the optics except the measurement of planets, we find that a single transit of Jupiter across the focal plane will measure the beam window functions to better than 0.3% for the channels at 100–217 GHz that are the most sensitive to the CMB. Constraining the beam with optical modeling can lead to much higher quality reconstruction. Conclusions. Depending on the optical modeling, the beam errors may be a significant contribution to the measurement systematics for n_s

Markov Chain Beam Randomization: a study of the impact of PLANCK beam measurement errors on cosmological parameter estimation

We introduce a new method to propagate uncertainties in the beam shapes used to measure the cosmic microwave background to cosmological parameters determined from those measurements. The method, which we call Markov Chain Beam Randomization, MCBR, randomly samples from a set of templates or functions that describe the beam uncertainties. The method is much faster than direct numerical integration over systematic `nuisance' parameters, and is not restricted to simple, idealized cases as is analytic marginalization. It does not assume the data are normally distributed, and does not require Gaussian priors on the specific systematic uncertainties. We show that MCBR properly accounts for and provides the marginalized errors of the parameters. The method can be generalized and used to propagate any systematic uncertainties for which a set of templates is available. We apply the method to the Planck satellite, and consider future experiments. Beam measurement errors should have a small effect on cosmological parameters as long as the beam fitting is performed after removal of 1/f noise.Comment: 17 pages, 23 figures, revised version with improved explanation of the MCBR and overall wording. Accepted for publication in Astronomy and Astrophysics (to appear in the Planck pre-launch special issue

The scalar perturbation spectral index n_s: WMAP sensitivity to unresolved point sources

Precision measurement of the scalar perturbation spectral index, n_s, from the Wilkinson Microwave Anisotropy Probe temperature angular power spectrum requires the subtraction of unresolved point source power. Here we reconsider this issue. First, we note a peculiarity in the WMAP temperature likelihood's response to the source correction: Cosmological parameters do not respond to increased source errors. An alternative and more direct method for treating this error term acts more sensibly, and also shifts n_s by ~0.3 sigma closer to unity. Second, we re-examine the source fit used to correct the power spectrum. This fit depends strongly on the galactic cut and the weighting of the map, indicating that either the source population or masking procedure is not isotropic. Jackknife tests appear inconsistent, causing us to assign large uncertainties to account for possible systematics. Third, we note that the WMAP team's spectrum was computed with two different weighting schemes: uniform weights transition to inverse noise variance weights at l = 500. The fit depends on such weighting schemes, so different corrections apply to each multipole range. For the Kp2 mask used in cosmological analysis, we prefer source corrections A = 0.012 +/- 0.005 muK^2 for uniform weighting and A = 0.015 +/- 0.005 muK^2 for N_obs weighting. Correcting WMAP's spectrum correspondingly, we compute cosmological parameters with our alternative likelihood, finding n_s = 0.970 +/- 0.017 and sigma_8 = 0.778 +/- 0.045 . This n_s is only 1.8 sigma from unity, compared to the ~2.6 sigma WMAP 3-year result. Finally, an anomalous feature in the source spectrum at l<200 remains, most strongly associated with W-band.Comment: 9 pages, 10 figures, 3 tables. Submitted to Ap

Reconstructing the shape of the correlation function

We develop an estimator for the correlation function which, in the ensemble average, returns the shape of the correlation function, even for signals that have significant correlations on the scale of the survey region. Our estimator is general and works in any number of dimensions. We develop versions of the estimator for both diffuse and discrete signals. As an application, we examine Monte Carlo simulations of X-ray background measurements. These include a realistic, spatially-inhomogeneous population of spurious detector events. We discuss applying the estimator to the averaging of correlation functions evaluated on several small fields, and to other cosmological applications.Comment: 10 pages, 5 figures, submitted to ApJS. Methods and results unchanged but text is expanded and significantly reordered in response to refere

The Q/U Imaging Experiment: Polarization Measurements of Radio Sources at 43 and 95 GHz

We present polarization measurements of extragalactic radio sources observed during the cosmic microwave background polarization survey of the Q/U Imaging Experiment (QUIET), operating at 43 GHz (Q-band) and 95 GHz (W-band). We examine sources selected at 20 GHz from the public, >40 mJy catalog of the Australia Telescope (AT20G) survey. There are ~480 such sources within QUIET's four low-foreground survey patches, including the nearby radio galaxies Centaurus A and Pictor A. The median error on our polarized flux density measurements is 30–40 mJy per Stokes parameter. At signal-to-noise ratio > 3 significance, we detect linear polarization for seven sources in Q-band and six in W-band; only 1.3 ± 1.1 detections per frequency band are expected by chance. For sources without a detection of polarized emission, we find that half of the sources have polarization amplitudes below 90 mJy (Q-band) and 106 mJy (W-band), at 95% confidence. Finally, we compare our polarization measurements to intensity and polarization measurements of the same sources from the literature. For the four sources with WMAP and Planck intensity measurements >1 Jy, the polarization fractions are above 1% in both QUIET bands. At high significance, we compute polarization fractions as much as 10%–20% for some sources, but the effects of source variability may cut that level in half for contemporaneous comparisons. Our results indicate that simple models—ones that scale a fixed polarization fraction with frequency—are inadequate to model the behavior of these sources and their contributions to polarization maps

Stacking catalog sources in WMAP data

We stack WMAP 7-year temperature data around extragalactic point sources, showing that the profiles are consistent with WMAP's beam models, in disagreement with the findings of Sawangwit & Shanks (2010). These results require that the source sample's selection is not biased by CMB fluctuations. We compare profiles from sources in the standard WMAP catalog, the WMAP catalog selected from a CMB-free combination of data, and the NVSS catalog, and quantify the agreement with fits to simple parametric beam models. We estimate the biases in source profiles due to alignments with positive CMB fluctuations, finding them roughly consistent with those biases found with the WMAP standard catalog. Addressing those biases, we find source spectral indices significantly steeper than those used by WMAP, with strong evidence for spectral steepening above 61 GHz. Such changes modify the power spectrum correction required for unresolved point sources, and tend to weaken somewhat the evidence for deviation from a Harrison-Zel'dovich primordial spectrum, but more analysis is required. Finally, we discuss implications for current CMB experiments.Comment: 10 pages, 7 figures, 2 tables, submitted to MNRA

Fast Pixel Space Convolution for CMB Surveys with Asymmetric Beams and Complex Scan Strategies: FEBeCoP

Precise measurement of the angular power spectrum of the Cosmic Microwave Background (CMB) temperature and polarization anisotropy can tightly constrain many cosmological models and parameters. However, accurate measurements can only be realized in practice provided all major systematic effects have been taken into account. Beam asymmetry, coupled with the scan strategy, is a major source of systematic error in scanning CMB experiments such as Planck, the focus of our current interest. We envision Monte Carlo methods to rigorously study and account for the systematic effect of beams in CMB analysis. Toward that goal, we have developed a fast pixel space convolution method that can simulate sky maps observed by a scanning instrument, taking into account real beam shapes and scan strategy. The essence is to pre-compute the "effective beams" using a computer code, "Fast Effective Beam Convolution in Pixel space" (FEBeCoP), that we have developed for the Planck mission. The code computes effective beams given the focal plane beam characteristics of the Planck instrument and the full history of actual satellite pointing, and performs very fast convolution of sky signals using the effective beams. In this paper, we describe the algorithm and the computational scheme that has been implemented. We also outline a few applications of the effective beams in the precision analysis of Planck data, for characterizing the CMB anisotropy and for detecting and measuring properties of point sources.Comment: 26 pages, 15 figures. New subsection on beam/PSF statistics, new and better figures, more explicit algebra for polarized beams, added explanatory text at many places following referees comments [Accepted for publication in ApJS

Probing Cosmology with Weak Lensing Minkowski Functionals

In this paper, we show that Minkowski Functionals (MFs) of weak gravitational lensing (WL) convergence maps contain significant non-Gaussian, cosmology-dependent information. To do this, we use a large suite of cosmological ray-tracing N-body simulations to create mock WL convergence maps, and study the cosmological information content of MFs derived from these maps. Our suite consists of 80 independent 512^3 N-body runs, covering seven different cosmologies, varying three cosmological parameters Omega_m, w, and sigma_8 one at a time, around a fiducial LambdaCDM model. In each cosmology, we use ray-tracing to create a thousand pseudo-independent 12 deg^2 convergence maps, and use these in a Monte Carlo procedure to estimate the joint confidence contours on the above three parameters. We include redshift tomography at three different source redshifts z_s=1, 1.5, 2, explore five different smoothing scales theta_G=1, 2, 3, 5, 10 arcmin, and explicitly compare and combine the MFs with the WL power spectrum. We find that the MFs capture a substantial amount of information from non-Gaussian features of convergence maps, i.e. beyond the power spectrum. The MFs are particularly well suited to break degeneracies and to constrain the dark energy equation of state parameter w (by a factor of ~ three better than from the power spectrum alone). The non-Gaussian information derives partly from the one-point function of the convergence (through V_0, the "area" MF), and partly through non-linear spatial information (through combining different smoothing scales for V_0, and through V_1 and V_2, the boundary length and genus MFs, respectively). In contrast to the power spectrum, the best constraints from the MFs are obtained only when multiple smoothing scales are combined.Comment: 19 pages, 9 figures, 5 table